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ABSTRACT 

Context. Our knowledge of circumstellar disks has traditionally been based on studies of dust. However, gas dominates the disk 
mass and its study is key to our understanding of accretion, outflows, and ultimately planet formation. The Spitzer Space Telescope 
provides access to gas emission lines in the mid-infrared, providing crucial new diagnostics of the physical conditions in accretion 
disks and outflows. 

Aims. We seek to identify gas emission lines in mid-infrared spectra of 64 pre-main-sequence stars in Taurus. Using line luminosities 
and other known star-disk-outflow parameters, we aim to identify correlations that will help to constrain gas heating, excitation 
mechanisms, and the line formation. 

Methods. We have based our study on Spitzer observations using the Infrared Spectrograph (IRS), mainly with the high-resolution 
modules. Line luminosities (or 3cr upper limits) have been obtained by fitting Gaussian profiles to the lines. We have further searched 
for correlations between the line luminosities and different parameters related to the star-disk system. 

Results. We have detected H 2 (17.03, 28.22 pm) emission in 6 objects, [Ne II] (12.81 pm) emission in 18 objects, and [Fe II] (17.93, 
25.99 pm) emission in 7 objects. [Ne II] detections are found primarily in Class II objects. The luminosity of the [Ne II] line (Lncii) is in 
general higher for objects known to drive jets than for those without known jets, but the two groups are not statistically distinguishable. 
L^eii is correlated with X-ray luminosity, but for Class II objects only. Ln c ii is also correlated with disk mass and accretion rate when 
the sample is divided into high and low accretors. Furthermore, we find correlations of L Ncn with mid-IR continuum luminosity and 
with luminosity of the [O I] (6300 A) line, the latter being an outflow tracer. L[ Fcn ] correlates with M ;lcc . No correlations were found 
between Lh, and several tested parameters. 

Conclusions. Our study reveals a general trend toward accretion-related phenomena as the origin of the gas emission lines. Shocks 
in jets and outflowing material are more likely to play a significant role than shocks in infalling material. The role of X-ray irradiation 
is less prominent but still present for [Ne II], in particular for Class II sources, while the lack of correlation between [Fe II] and [Ne 
II] points toward different emitting mechanisms. 

Key words. ISM: jets and outflows - Infrared: stars - Protoplanetary disks - Stars: formation - Stars: pre-main sequence - Stars: 
protostars 

>^ ' 1. Introduction protoplanetary disks, namely thermal emission from dust grains 

H \ in the millimeter and the infrared, and molecular line emission 

W , Protoplanetary disks are formed by a mixture of gas, dust and from co in the su b-millimeter, typically disappear after a cou - 
ices. The bulk of their mass is in gaseous form (mainly H 2 ) D l e of Mvr (e.g JDuvert et alJ200o[lAndrews & Williamsl2005h . 
but the opacity is dominated by dust at infrared (IR) and opti- However, it is not clear whether the material has been completely 
cal wavelengths. Ro-vibrational transitions of H 2 in the IR are c i ea red out or if present observations are not sensitive enough to 
difficult to stimulate due to the combination of low tempera- allow the detection of its remnant. Recent evidenc e indicates that 
tures found in protoplanetary disks and the lack of a permanent the gas and dust lifetime s are not the same (e.g., iPascucci et al] 
electric dipole moment of the H 2 molecule. For this reason, our |2009; Fed ele et al.ll2oToh . Detections of quiescent emission of 
current knowledge of protoplanetary disks is largely based on mo i ecu iar H 2 in weak-lined T Tauri stars (WTTS) that do not 
the study of dust (e.g Beckwith et al.lll990t |Hartmann||2008t show i n f rare d excess, i. e., those that are not expected to have 
| Mann & Williams ||2009|). a disk, sup port the hypothesis that gas lives longer than dust in 

One of the key questions in star formation is the disk life- disks (e.g.. lBarv et al.ll2002ll2008ah . To address the question of 
time, since it provides a strong constraint on the planet forma- disk dissipation timescales, studies of the gas component of pro- 
tion time scale. Many stud ies have been dedicated to this prob- toplanetary disks independent of dust are particularly needed, 
lem (e.g., Carpenter et al. 20091 lMevetlf2 009). The signatures of 



1 



Baldovin-Saavedra et al.: Searching for gas lines in the young population of the TMC 



To study gas in protoplanetary disks, molecular spec- 
troscopy, in particular emission lines, can be used. As the tem- 
perature and density of the gas vary as function of the radius 
in the disk and the mid-plane distance, spectral lines at differ- 
ent wavelengths probe different disk regions, and gas at differ- 
ent physical conditions. In the sub-millimeter and millimeter 
regimes, observations probe cold gas (T ~ 20 - 50 K) at large 
radii (> 20 AU) through emission lines such us CO speci es and 
HCO + (e.g jMannings & Sargentll997tlPietu etai]|2007l) . 

In the near-IR, the H2 emission arises fro m ro-vibrational 
transi tions, probing gas at T~ 1000 K (e.g ., [Weintraub et al.l 
2000; ICarmona et alj|2007t IBarv et al.ll2008bl) . After H 2 , CO is 
the most abundant species found in disks. The different excita- 
tion conditions of this molecule probe a wide range of disk radii; 
the first overtone band at 2.3 pm traces hot (~ 4000 K) and dense 
material (n > 10 10 crrT 3 ) of the innermost part of protoplanetary 
disks. Observations show that the emissio n has the character- 
istic shape suggestive of a rotating disk dCarr| [l989:). The CO 
ro-vibrational band at 4.7 Ami is estimated to arise within 1-10 
AU from the star (e.g..lNaiita et al.l2003l;lBlake & Boogertl2004t 
iRetfig et al.ll2004tlSalvk et al.ll2007l;ICarmona et al.ll2008bl) wife 
characteristic temperatures between 100 to a few 1000 K. 

The mid-IR traces the warm gas at radii of a few AU 
up to several tens of AU (the giant planet formation re- 
gion). At these wavelengths, H2 emission is observed, aris- 
ing from pure rotational transitions at 12.28, 17.03, and 
28.22 /mi (e.g.. Ivan den Ancker et aTlll999h ILahuis et alj|2007l 
iMartin-Zaidi et al.ll2007LlBitner et al]l2008l) . Furthermore, many 
detections of forbidden atomic lines have been reported, such as 
[Ne II] (12.81 /mi), [Fe I] (24.0 /mi), [Fe II] (17.93 and 25.99 
/mi), and [S III] (18.71 /mi), opening up a new window into the 
diagnosis of gas in protoplanetary disks. 

When studying the gas emission from protoplanetary disks, 
it is important to discriminate among emission produced by (i) 
gas in jets or outflows, (ii) accreting gas falling onto the star, and 
(iii) gas in the disk itself. High spectral and spatial resolution 
observations enable the study of the line profiles and possible 
velocity shifts with respect to the stellar rest velocity. However, 
the heating mechanism of g as and the ionization of atomic 
species is not well understood. Glassg oTd et alJ J2007) predicted 
that stellar X-rays heat the upper layers of protoplanetary 
disks within 25 AU from the star, exciting atomic species. 
Gorti & Hollenbachl (120081) considered both stellar X-rays 
and ultraviolet (UV) irradiation, proposing several emission 
lines as possible diagnostics of the gas in disks irradiated by 
high energy photons: forbidden lines of [Ne II] (12.81 /mi), 
[Ne III] (15.55 /mi), [Fe II] (25.9 9 urn), [O I] (63 u rn), and pure 
rotational lines of H 2 . Moreover. lAlexa nder (2008) successfully 
reproduced the observed [Ne II] line luminosities by modeling 
the line profile using a photoevaporative wind model, which 
considered the UV flux from the central star as the only 
ionizing source. In addition, far-uv (FUV), extreme-UV (EUV), 
and X-ray photons may be produced by accretion shocks as 
the m aterial falls into the central star dHollenbach & GortH 
2009). Forbidden emission lines in the optical and near-IR 
have been used for a l ong time to probe the shocked gas 
from T Tauri stars (e.£ 

20091: lAgra-Amboage et al.l 120091) . 
studied the infrared emission from 



Hartigan & Hillenbran 
Hollenbach eta!T "(l'l989 



interstellar shocks and proposed that for low pre-shock densities 



(n < IO 4 - 3 cnT 3 ), forbidden emission lines of [Fe II], [Ne II], 
[Si II], and [O I] are particularly strong. 

Forbidden emission lines are used to trace the electron den- 
sities of the environment where they originate; typical densities 



of 10 3 - 10 7 crrT 3 are probed. At these densities, the origin is in 
low-density winds or shocks rather than in the infall of material 
through accretion funnels. On the other hand, permitted emis- 
sion lines exhibiting large widths an d inverse P-Cygni fea tures 
are associated with infalling material dMuzerolle et al.l l2001). 

Up to now, observations have not been conclusive in giv- 
ing a broadly accepted explanation of the heating and ioniz- 
ing mechanism exciting the gas component in protoplanetary 
disks. We aim to contribute to the understanding of this prob- 
lem by studying a large sample of pre-main sequence stars of 
low and intermediate mass. We have used observat ions done 
by the Spitzer Spac e Telescope ( Werner et alJ 120041) Infrared 
Spectrograph (IRS, iHouck et al.l 120041) in order to study the 
emission from gas lines. Although higher spectral resolution ob- 
servations are needed in order to obtain kinematical information 
for the emitting gas, the high sensitivity of Spitzer does enable us 
to access sources not detectable from the ground. Furthermore, 
the high number of sources observed allows the extraction of sta- 
tistically significant information from a large sample. We have 
made use of archival spectra of 6 1 pre-main sequence stars plus 
3 objects from our own Spitzer program. This study is comple- 
mentary to previous investigations; our work includes objects of 
different evolu tionary stages, w hile the somewhat related previ- 
ous study by Gtide l et alJ (1201 Oh concentrated on a pure Class II 
sample. Additionally, in our analysis we have made a distinction 
based on the YSO class, in contrast to previous studies (e.g., 
ILahuis et al.ll2007l:lFlaccomio et al.ll2009l) . 

Sections|2]and[3]are dedicated to the description of the sam- 
ple and of the data analysis, respectively. Results are presented 
in Section [4] and discussed in Section [5] together with the pre- 
sentation of correlation tests. A brief summary of the results and 
conclusions are given in Section [6] In Appendix lAl (available in 
the electronic version) we have included additional information 
on individual objects whenever one of the studied gas lines was 
detected. 



2. Sample 

Our study is based on objects from the Taurus Molecular Cloud 
(TMC). This star-forming region is known to form mainly low- 
mass stars, allowing the study of the star formation process with- 
out the influence of high-mass young stars, which are known 
to have strong UV radiation fields. Its nearby distance (137 pc, 
iTorres et"ai1l2007l) enables us to access fainter sources. All the 
objects have similar ages, lie at nearly the same distance, and 
are presumably formed out of the same parent cloud. 

A further motivation to focus our study on the TMC re- 
gion is the large amount of multi-wavelength observations avail- 
able. The TMC was indeed been observed extensively in X- 
rays with XMM-Newto n: the XEST surve y covered ~5 square 
degrees of the TMC, iGiidel et all 120071 In the optical, the 
Canada-France-Hawaii Telescope (CFHT) covered 28 square 
degrees of the TMC (Monin et al. 2010, in preparation), and 
the Sl oan Digital Sky Survey (SPSS) covered 48 square de- 
grees (Padmanabhan et al. 2008). In the IR, Spitze r mapped ~43 
squar e degrees (Spitzer Taurus Legacy Project; Padge ttet al.l 
2007, 2011, in preparation). Multiwavelength observations are 
crucial in order to test the influence of the stellar X-rays in the 
heating and excitation of the gas species. 

Our study is based on publicly available archived data ob- 
tained by IRS with its high resolution modules (short-high - 
SH - and long-high - LH), with a resolving power of R~ 600 
(~ 500 km s _1 ) for wavelength coverage between 9 and 37 /mi. 
Most of the targets studied belong to GTO (Guaranteed Time 
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Observations) programs; several publications exist that focussed 
on th e study of the silicate and ice features (e.g., IWatson et alJ 
12004 iFurlanet al.1l2006t iFurlan et al.1l2008l) . In particular, we 
used data from the GTO program number 2, from the Legacy 
programs 171, 179, and from the GO (General Observer) pro- 
grams 3470, 3716, and 30300. In addition, we included 3 low- 
mass pre-main sequence stars identified by the Spitzer TMC sur- 
vey for which we obtained IRS high and low-resolution spec- 
tra (program 50807, PI. Baldovin-Saavedra). For one object 
(SST 041412.2+280837) we only have the spectrum taken with 
the low-resolution modules (long-low - LL - and short-low - 
SL). 

Our sample consists of 64 low and intermediate mass pre- 
main sequence stars. Their infrared classification is taken from 
the literature and is based on the slope of the spectral energy dis- 
tribution (SED) in the infrared. The sample includes 10 Class I 
objects, 40 Class II, 4 Class III, and 6 intermediate mass Herbig 
AeBe stars. We found 4 objects classified in the literature as both 
Class I and Class II. Table Q] (available in the electronic version) 
summarizes the properties of the stars in our sample, including 
the bibliographic references for the adopted values. 

3. Data Analysis 

3.1. Initial Steps 

The data were retrieved from the Spitzer archive; we started 
the analysis from the post-basic calibrated data, i.e., spectra ex- 
tracted automatically by the pipeline (version S 15.3.0 for SH and 
S17.2 for LHQ). The spectra of the three targets from our own 
program were treated using the pipeline version S18.7.O0 

In a first step we used IRSCLEAN0 to create masks for the 
rogue pixels from the basic calibrated data images (rogue pix- 
els have abnormally high dark current and/or photon respon- 
sivity). Once the images were free from obvious rogue pixels, 
we re-extracted the spectra using SPICE (Spitzer IRS Custom 
Extraction) regular extraction. Furthermore, in the cases where 
we detected emission lines or where the detections were not ob- 
vious, we also used optimal extraction (details on the regular 
and optimal extraction can be found in the SPICE User's Guide 
Version, 16.2. 1+). 

Some of the spectra also presented residual fringes; these 
interference patterns are typical for infrared detectors as their 
surface coating acts as a Fabry-Perot etalon, leading to a sinu- 
soidal modification of the detector response of spectrally dis- 
persed light. Although they were corrected during the pipeline 
processing, residual fringes remained in a few cases that we re- 
moved using the tool IRSFRINGE 3 . 

3.2. Background 

The mid-IR background comes predominantly from zodiacal 
light at the IRS wavelengths, while diffuse emission from the 
Taurus cloud contributes only little. Instrumental background 
(rogue pixels) is also significant for the high-resolution mod- 
ules. Therefore, sky observations are needed in order to correct 
for the sky and instrumental background. In this study, however, 

1 We checked for a few sources that the latest version of the pipeline 
(SI 8.7.0) did not change the obtained line luminosities. 

2 Our observations were reduced from the basic calibrated data, cor- 
recting for the co-adding bug in SI 8.7.0 and for time and position- 
dependent dark currents on the LH array (DarkSettle 1.1). 

3 Available from http://ssc.spitzer.caltech.edu 



most of our source sample had no dedicated background obser- 
vations. Good rogue pixel removal proved to be crucial for our 
emission line study. 

Some of the GTO observations were performed using a se- 
ries of 6 pointings per object with a shift corresponding to ap- 
proximately half of the PSF In these cases, for the on-source 
observations we used the 2 pointings with higher flux and aver- 
aged them; these pointings are equivalent to the standard staring 
mode observations, i.e., 2 slit positions with the target located at 
30 and 60% of the slit length. For each object, we checked us- 
ing Spitzer Infrared Array Camera (IRAC) 8 pm images that the 
pointings correspond to the SIMBAD position of the target, and 
that the flux levels of the 2 used pointings were equivalent. 

One central issue lies in understanding the origin of the emis- 
sion. It could be that the lines detected are not really emitted by 
the star-disk system, but rather come from the background or 
a nearby source. Our data are in the form of an echellogram, 
which is a two-dimensional representation of a spectrum con- 
taining both spatial and spectral information. As a first step, we 
searched each spectrum and searched for evident line emission; 
when found, we examined in the echellogram two spectral re- 
gions with no emission lines adjacent to the position of the line 
studied. These two regions were used as a measure of local back- 
ground plus stellar continuum; they were averaged and then sub- 
tracted from the line region. In this way, any possible extended 
emission can be identified; moreover, by removing the contin- 
uum, the detection becomes more evident. This method did not 
result in the discarding of any detections. However, it is not sen- 
sitive to potential extended emission from the direct vicinity of 
the star. 

3.3. Line Fitting 

In order to obtain line luminosities, we used a function consist- 
ing of a Gaussian plus a linear component to account for the 
continuum within a range of typically +0.6 pm of the line. All 
the parameters of the function (height, center, and width) were 
left free to vary. The integrated flux was calculated as the area 
under the Gaussian fit once the continuum was subtracted, and 
the line luminos ity was calculated u sing a distance to the TMC 
of 137 pc from iTorres etall (l2007lP l A detection was claimed 
when the peak of the line was higher than 3 times the level of 
the noise in the spectra. The noise level is the standard deviation 
of the continuum flux (cr) in the wavelength range of interest. 

Most of the detected lines were spectrally unresolved; the 
full width at half-max (FWHM) obtained from the fit was similar 
to the instrumental FWHM (~ 500 km s _1 ). In the cases where 
the lines were not been detected, we obtained 3<x upper limits 
for their luminosities by multiplying the 3cr noise value of the 
continuum by the FWHM. 

4. Results 

Our study is focused on the search for pure rotational emission 
lines from H2 (12.28, 17.03, 28.22 pm), and forbidden emis- 
sion lines from [Ne II] (12.81 pm), [Ne III] (15.56 pm), [Fe II] 
(17.94, 25.99 pm), and [S III] (18.71, 33.48 pm). We detected 
[Ne II] in 18 objects in our sample, [Fe II] in 7 objects, and H2 
in 6 objects. We did not detect [Ne III] or [S III]. The line lu- 
minosities and upper limits are presented in Table [2j with the 

4 The distance to the TMC has often been considered to be 140 pc. 
In comparison to the 137 pc used in this study, this can result in a dif- 
ference of the order of 4% in the line luminosities. 
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luminosities expressed in units of 10 28 erg s _1 ; bold fonts in- 
dicate detected lines, and normal fonts are 3<x upper limits. In 
Appendix lAl (available in the electronic version) we provide in- 
formation on the single objects with detected emission lines. 

In Fig. Q] we present the spectra of the targets with detected 
[Ne II] line emission, whereas we present the spectra surround- 
ing the position of the [Ne II] line for the sources with non- 
detections in Appendix 151 (available in the electronic version). 
Figures [2] and [3] show the line spectra of the sources with de- 
tected emission lines from H2 and [Fe II]. In all figures, the red 
line shows the Gaussian fit to the line, the dashed line shows 3<x 
detection threshold, and the vertical dashed-dotted line shows 
the expected position of the [Ne II] line. The title of each panel 
is the name of the target. 



5. Discussion 

5.1. [Ne II] emission 

The search for [Ne II] emission has gained importance in the last 
few years thanks to an increasing number of reported detections 
in the spectra of young stars observed ma i nly with Spitzer (e.g. , 



Pascucci et al. 2007; 



Lahuis et al 



Flaccomio et al.ll2009l:lGudel et ali bOlO). The sample presented 



20071: lEspaillat et al] 120071: 



in Pascucci et al. (2007) included six targets showing evidence 
of optically thick di sks. The [Ne II ] line was detected in four 
of them. Moreover, iPascucci et al.1 d2007l) reported a correla- 
tion between the lu minosity of the [Ne II] line and the stel- 
lar X-ray luminosity. lLahuis et alj d2007l) detected [Ne II] in 15 
out of 76 low-mass pre-main sequence stars from different star 
forming regions, and also reported one detection of [Ne III]. 
Those authors concluded that the [Ne II] detections and the 
[Ne II]/[Ne III] line flux ratio are consistent with excitation 
through X-ray heating of the circumstellar disk from the central 
star, althou gh excitation through E UV might contribute as well. 
In addition, Fl accomio et al.l {2009) obtained Spitzer IRS spectra 
for 28 young stellar objects in the Ophiuchus region, obtaining 
10 detections of [Ne II] and one detection of [Ne III]. They did 
not find a correlation between the luminosity of the [Ne II] line 
and the stellar X-ray luminosity. Instead, the luminosity of the 
[Ne II] line seemed to be correlated with disk mass accretion 
rate, consistent with lEspaillat et"ai1 (120071) . iGiidel et al.1 <l2010h 
based their study on 92 sources observed with Spitzer IRS, ob- 
taining detections of [Ne II] in 58 of them. A weak correlation 
was found between the luminosity of the [Ne II] line and the stel- 
lar X-ray luminosity, whereas a good correlation was also found 
between L[N e ii] and jet/outflow parameters, such as the luminos- 
ity of the [O I] line. 

We have detected the [Ne II] line in 18 objects, represent- 
ing 28% of our sample. This de tection rate is similar to the 
results of Flaccomi o et al ] (120091) . That study was based only 
in the Ophiuchus regio n, obtaining a det ection rate of 36%. 
lLahuis et al.l ((2007) and Gii del et al.l d2010l) used a mixed sam- 
ple with stars from different star-forming regions, obtaining 
~ 20 and ~ 60 % detections respectively. The sample from 
IGiidel et al.l J20T0I) did not contain any Class I or Herbig AeBe 
stars. Furthermore, their data red uction differed from ours, as 
they used the method introduced in lLahuis et al.l (12010 ) in which 
the compact and extended component of the line emission are 
separated during the spectral extraction procedure. 

Among the objects in our sample with [Ne II] detections, 
the distribution in terms of evolutionary state is the following: 
1 1% of the [Ne II] detections correspond to Class I objects, 67% 
of the detections correspond to pure Class II, and 22% of the 



detections correspond to objects classified as Class I and II in 
the literature. 

The highest luminosity of the [Ne II] line in our study came 
from T Tau, which is also one of the most luminous young 
star in the X-rays (Lx — 8 x 10 30 erg s _1 ). This system is 
know n to drive an outflow responsible for the [Ne II] emis- 
sion (IVanBoekel et alj 12009). L1551 IRS 5 has also a [Ne II] 
luminosity comparable to that of T Tau, of order 10 30 erg s . 
Previous studies in the optical and near-IR have shown that 
LI 551 IRS 5 drives a strong molecular outflow (see Appendix 
lAl . so there is a suggestion in this case that the [Ne II] could 
be produced in an outflow/jet environment. We note that its X- 
ray luminosity is not well constr ained, given that the source is 
very absorbed dBallv et al.ll2003h : the value available is associ- 
ated with the jet. For this reason, in the correlation tests pre- 
sented in Section l5.1.2l we do not consider the X-ray luminosity 
ofL1551IRS5. 

Other sources with [Ne II] detections known to drive out- 
flows are AA Tau, DG Tau, DG Tau B, UY Aur, and CoKu Tau 1 . 
In general, these sources show high [Ne II] luminosities (> 
10 29 erg s- 1 ) except for AA Tau (L [Ne n] ~ 3 x 10 28 erg s" 1 ). We 
might conclude that for high [Ne II] luminosities, the heating of 
gas and excitation of neon is caused by the outflow activ i ty. Thi s 
trend is actually observed in the study of Gii del et al.l yOlO); 
sources known to drive jets or outflows show an excess of [Ne II] 
emission when compared with sources that do not present jets. 

To check the suggested difference between sources driving 
jets and sources without jets, we have divided our sample into 
sources showing jet emission (typically forbidden line emis- 
sion in the optical and near-IR) and sources without known jets. 
Using the Kaplan- Meier (K-M) estimator available within th e 
ASURV0package dFeigelson & Nelson! 1985tlisobe et al.ll986h . 
we have extracted the cumulative distribution of [Ne II] lumi- 
nosities for these two groups. The K-M estimator is the ba- 
sis of survival analysis. It is a generalized maximum-likelihood 
estimator for the population from which the sample is drawn. 
Although it was first used in right censoring problems (i.e., lower 
limits), it has been adapted for application to left censoring prob- 
lems (i.e., upper limits). The basic concept underlying the K-M 
estimator is the survivor function, which is defined as the prob- 
ability of a value being larger than or smaller than the upper or 
lower limit, depending on the type of censoring. The survivor 
function is a non increasing step function which only jumps at 
uncensored values. When the sample is large, it approaches the 
true survivor function of th e population. For details we refer to 
Feige lson & Nelson! (fl98l . 

The sub-sample of sources without jets has 38 objects; 
among them, 9 have [Ne II] detections. The mean of the dis- 
tribution is L[Ne ii] — (0.8^03) x 10 28 erg s -1 . The sub-sample 
of sources with jets has 26 objects; among them we have 9 
detections of the [Ne II] line. The cumulative distribution of 
[Ne II] luminosities is extended towards higher values, with a 
mean of L[N e ii] = (1.6*9 - |) x 10 28 erg s _I , almost twice the 
mean value of the sub-sample of sources with no jets. However, 
both distributions overlap. Moreover, the median of the two sub- 
samples are very close, with L[N e ii] = 4 x 10 28 erg s _1 and 
L[Ne ii] = 5 x 10 28 erg s _1 for jet and non-jet sources, respec- 
tively. Thus, w e cannot distinguish between the two populations, 
in contrast with lGiidel et al.l d2010). Th is difference might be due 
to the sample selection criteria. In fact. IGiidel et al.l (l2010) chose 
a sample composed only of Class II objects. Furthermore, their 
sub-sample of jet-driving sources, unlike ours, did not contain 



5 Available from http://astrostatistics.psu.edu/statcodes/asurv 
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Fig. 1. Detections of [Ne II] at 12.81 /urn. The spectra are plotted as solid histogram, the lines were fitted using a gaussian profile 
plus a polynomial component for the continuum, the fit is plotted using a solid line. The 3<x detection threshold is plotted as a dashed 
line. The name of each star is labeled in the top of each plot. We recall that the spectrum of SST 041412.2+280837 was obtained 
using the low resolution modules (SL and LL). 



objects with indirect evidence of jets (i.e., high-velocity emis- 
sion line components). 

The number of detections of [Ne II] in our sample (18 ob- 
jects) has prompted us to search for possible correlations or 
trends between the luminosity of the line and parameters related 
to the star-disk system, and outflow/jet indicators. The parame- 
ters studied are: stellar mass (M*), stellar luminosity (L star ), X- 
ray luminosity (Lx), Ha equivalent width, mass accretion rate 



(M acc ), wind mass loss rate (M w j n d), luminosity of the [O I] line 
in the optical (Ljoij), and luminosity of the [Fe II] line at 26 jum 
(L[Feii])- If correlations exist, they would give new hints on the 
neon excitation mechanism. 

Due to the high number of upper limits for the luminosity of 
the [Ne II] line in our sample, we have tested the correlations us- 
ing survival analysis implemented in the ASURV package. We 
have used 3 statistical tests: Cox's proportional hazard model, 
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Fig. 2. The plotted spectra show the regions around the H2 lines at 12.28, 17.03 and 28.22 fim for the sources where we have detected 
at least one of the lines. 



generalized Kendall's r rank correlation, and Spearman's p. The 
Kendall and Spearman methods allow the use of upper limits in 
one or both of the variables tested. Therefore, we have included 
upper limits for all the parameters tested. In contrast, the Cox 
hazard model allows only upper limits in the dependent vari- 
able. The Spearman method is not reliable for data sets smaller 
than 30 data points. Nevertheless, we report for completeness 
the probabilities obtained with this method even when the num- 
ber of data points is smaller than 30. Calculations were done 
using the logarithm of values. Table [1] summarizes our results, 



in which we have included the values for the correlation param- 
eters between parentheses for completeness. In the case of the 
Cox hazard method, the parameter reported is the x 2 of the re- 
gression fit. For the Kendall method, the value reported is the 
ratio between the correlation parameter r and its standard de- 
viation <j T , a function of the correlation. For the Spearman's p 
method, the parameter reported is the estimated correlation. 

We have tested the correlations using the whole sample (i.e., 
all infrared classes). In addition, we have created a sub-sample 
with the Class II sources, given that the majority of the [Ne II] 
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detections are from these objects. In Figs. [4] and E] we present 
L[Ne ii] as a function of the different parameters studied for the 
complete sample and for the sub-sample of Class II objects sepa- 
rately. Although we find significant correlations (at a 5% thresh- 
old) based on several methods included in the ASURV package, 
there is a certain chance that the conclusion of correlation is in- 
correct. To this end, we have used a bootstrap method to de- 
termine the chance of finding a correlation by chance. Such a 
method implies the creation of a sampling distribution from the 
original data. In practice, we created 5 x 10 4 random samples in 
each testing case by keeping the independent variable unchanged 
and permuting each time the dependent variable (L[Ne ii] in this 
case) in order to break the links between each pair of points, 
keeping the number of detections and upper limits unchanged. 
We have applied the correlation tests using the ASURV package 
for each sample. Then we have constructed the distribution of the 
correlation parameters (i.e., t for Kendall, p for Spearman, and 
X 2 for Cox), and we have looked into the percentage of cases 
with correlation parameters larger than the ones obtained for 
the original sample. In all correlations tested, the probabilities 
were generally less than 4%, thus confirming the correlations 
presented in Table Q] 

5.1 .1 . Correlations with stellar properties 

We have tested correlati ons with the stel l ar mas s and luminos- 
ity. In a previous study, Flaccomi o et alJ {2009) found that the 
stellar mass does not influence the [Ne II] luminosity. In our 
sample, the detected [Ne II] lines are concentrated in sources 
at the low end of the mass range (with stellar masses typically 
<1 Mp) ; only two object s with M* > 1 M show [Ne II] emis- 
sion. In Flaccomio et al. (2009), the detected [Ne II] lines are 
spread across a larger range of stellar masses. Despite the dif- 
ferent mass ranges probed by the two samples, the result is the 
same: L[N e ii] is not dependent on the stellar mass. Since the stel- 
lar mass and lu minosity are related through the relation L* oc Mj 
dHartmannll2009h . we expect to obtain a similar lack of correla- 
tion between L[N e ii] and the stellar luminosity, which is indeed 
confirmed. 



5.1.2. Correlations with high-energy irradiation 

The stellar X-ray luminosity is one of the key parameters driv- 
ing the [N e II] emiss i on pro posed in several theoretical studies, 
e.g., [Glassgold et al] d2007l) . T he first observational study sug- 
gesting such a correlation was iPascucci et al.l {2007), although 
this result was based on a sample of just fo ur objects, spanning 
a small range in Lx. Espaillatet al. (2007) increased the sam- 
ple to seven objects but did not find any evident correlation be- 
tween LrN e m and Lx. T he same result was obtained recently by 
iFlaccomio et al 1 (120091) based on a larger sample and spanning a 
larger range in Lx- Nevertheless, IGudel et all 1201 Ch have found 
a statistically significant dependence on the stellar Lx for their 
sample confined to Class II objects and spanning a broad range 
ofL x . 

Our sample spans a range of X-ray luminositis between 
10 29 and ~ 10 31 erg s _1 , similar to the sample presented by 
Giidel et al.l (12010). The correlation between LNeii and Lx is non- 
existent for the complete sample (a ll infrared classes), i.e., con- 
sistent with Flaccomio et ID {2009). But we observe a correla- 
tion with Lx for the su b-sample of Class II objects, consistent 
with Giidel et al] d2010t) . These results are not contradictory: in 
fact, the difference between the previous studies lies in the sam- 



ple selection. IFlaccomio et al ] (120091) used a mixed sample, in- 
cluding Class I to Class III objects, similar to our sample (al- 
though we have included also Her big AeBe stars), whereas the 
sample used in IGudel etaD d2010j) contains only Class II ob- 
jects. Therefore, we do expect to find a lack of correlation like in 
Flaccomio et a II (120091) when usin g our whole sample, and the 
same result as in IGudel et al.l (12010) when restricting the sample 
to Class II objects. 

5.1.3. Correlations with accretion parameters 

The mass accretion rate (M acc ) is an indirect measurement of the 
evolutionary state of the sources. We expect M fc to decrease 
as the star evolves toward the m ain sequence. Espaillat et al.l 
(2007) and Fla ccomio et al (l2009l) found a posi t ive co rrelation 
between L[N e nj and M acc . Hollenbac h & GortH (l2009l) studied 
the irradiation of the disk by high-energy photons (EUV and 
soft X-rays). They showed that in order for the high-energy pho- 
tons to penetrate the strong protostellar winds and illuminate 
the disk beyond 1 AU, low wind mass-loss rates were needed 
(M w j n d < 10~ 9 Mpyr" 1 ). T his critical value was translated in 
iHollenbach & Gortil (120091) into a critical disk mass accretion 
rate of M acc < 10~ 8 M yr~'. Sources with low mass accretion 
rates would be more suitable to study the effects of high-energy 
photons. Following this idea, we have split our sample between 
high and low acc retors, setting the lim i t at M acc = 10~ 8 M yr _1 , 
as proposed in Hollen bach & Go rti (2009). We have, there- 
fore, defined two sub-samples: M acc _hi e h ^ 10~ 8 M Q yr _1 , and 
M acc _i ow < 10~ 8 M Q yr _1 . For the sub-sample of Class II ob- 
jects, we find a positive correlation between L[N e ii] and M acc 
for the group M acc _hi g h and a negative correlation with L[N e ii] 
for the group M acc _i ow . This result is not found when consider- 
ing the sub-sample of Class II objects with no separation in the 
mass accretion rate range. We observe a difference between low 
and high accretors; the anti-correlation found for low-accretors 
is statistically more significant. But we are cautious about this 
result: this correlation has been discussed in IGudel et al.l (l2010) 
in the sense that jet driving sources are younger and, in general, 
have higher mass accretion rates than sources without jets. They 
argued that the segregation between both samples might produce 
an apparent correlation. Indeed, in our study, the sub-sample 
of low-accretors (M acc _i ow ) is dominated by non-jet sources; 18 
out of the 21 sources in this sub-sample do not have evidence 
of jets. On the other hand, in the sub-sample of high accretors 
(M acc _hi g hX which accounts for 27 sources, 16 present observa- 
tional evidence of jets/outflows. 

The Ha line (seen usually in emission in young stars) is 
another tracer of accretion. Its equivalent width is commonly 
used to differentiate between accreting and non-accreting stars. 
Typically, EW(Hq') > 10 A for accreting stars, and EW(Ha)< 
10 A for non-accreting stars (lMartinlll997l) . The result for the 
test between L[N e nj and EW(Ho-) is not clear in our sample. 
Although we found a correlation for the whole sample (i.e., all 
infrared classes), this result is not repeated for the sub-sample of 
Class II objects. 

Finally, we have tested a correlation of L[N e ii] with the lumi- 
nosity in the mid-IR, in order to check if bright emission lines are 
related to a bright continuum emission. The mid-IR luminosity 
has been obtained for each object by integrating the continuum 
flux in the SH range (9 to 19 /urn). Since in most cases we do not 
have background observations, we have over-plotted the SH and 
SL spectra and checked that the continuum emission is domi- 
nated by the source and not by the background. Given the small 
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Table 1. Results from survival analysis correlation tests for the luminosity of the [Ne II] line. 



All infrared classes 

N° Upper Limits Probability and Correlation Parameter 



Variable 


N° Points 


Total 


[Nell] 


Cox 


Kendall 


Spearman 


M, 


53 


38 


38 


0.92 (0.010) 


0.62 (0.49) 


0.55 (0.084) 


L, 


57 


42 


42 


0.56 (0.34) 


0.82 (0.23) 


0.93 (0.012) 


Lx 


47 


35 


33 




0.09 (1.7) 


0.18 (0.18) 


Lir 


64 


46 


46 


0.04 (4.4) 


0.01 (2.4) 


0.01 (0.31) 




56 


42 


42 


0.04 (4.3) 


0.03(2.1) 


0.01 (0.3) 


M acc 


47 


36 


32 




0.62 (0.50) 


0.37 (0.13) 


M acc _high 


27 


21 


20 




0.13 (1.5) 


0.04 (0.39) 


Macc-low 


19 


15 


12 




0.20 (-1.3) 


0.55 (-0.14) 


L[on 


27 


21 


16 




0.02 (2.2) 


0.02 (0.46) 


LxL[oi] 


22 


17 


12 




0.06(1.9) 


0.03 (0.46) 


LxM w j n d 


36 


28 


23 




0.67 (0.42) 


1.00 (0.00) 


L[FeII] 


64 


61 


46 




0.09(1.7) 


0.10(0.21) 


Class II 


M, 


40 


26 


26 


0.10(2.6) 


0.46 (-0.74) 


0.27 (-0.17) 




40 


26 


26 


0.94 (0.005) 


0.85 (0.19) 


0.88 (0.02) 


Lx 


35 


24 


22 




0.04 (2.0) 


0.05 (0.34) 


Lir 


44 


28 


28 


0.03 (4.7) 


0.004 (2.8) 


0.01 (0.39) 


EW(Ha) 


41 


28 


28 


0.26(1.3) 


0.29(1.1) 


0.23 (0.19) 


M acc 


36 


26 


24 




0.90 (0.12) 


0.99 (-0.002) 


M acc _high 


19 


14 


13 




0.05 (1.9) 


0.04 (0.49) 


^Lcc-low 


17 


12 


10 




0.02 (-2.2) 


0.05 (-0.49) 


L[ rj 


24 


18 


18 




0.02 (2.3) 


0.03 (0.45) 


LxL[oi] 


19 


14 


13 




0.10(1.9) 


0.09 (0.44) 


LxM w i„d 


30 


22 


19 




0.61 (-0.51) 


0.63 (-0.90) 


L[FeII] 


44 


42 


28 




0.62 (0.50) 


0.41 (0.12) 



Notes. The first column shows the variable for which the correlation is tested, the second column shows the total number of points considered, the 
third column shows the total number of upper limits used (dependent and independent variable), the fourth column shows the number of upper 
limits for the luminosity of the [Ne II] line (dependent variable), and the probabilities (P) of a correlation achieved by chance are presented in 
columns 5, 6, and 7, where P < 5% are shown in bold fonts. The value of the correlation parameters are shown in parentheses. We have defined 
M acc _hi g h and Macc-iow as follows: M acc _ hii;h > 10~ 8 M o yr~' and M ilcc -i ow < 10~ 8 M Q yr~'. The Cox method cannot be used when upper limits are 
present in both of the tested variables. See text for more details. 



aperture and wavelength coverage of the SH modu le, the back- 
groun d contribution is expected to be small (e.g., IFurlan et alJ 
2006). We have obtained a correlation between L[N e ii] and Lir. 
This result is confirmed by the bootstrap technique. The origin 
of this correlation is not clear, although not likely to be the prod- 
uct of bias in the observations. The rms of the continuum in- 
creases with Lir, but the signal-to-noise ratio for the detected 
[Ne II] lines is distributed uniformly over the Lir spanned by 
our sample, with no preference for either high or low values 
of the signal-to noi se ratio. This correlation was also found by 
lLahuis et alJ d2010h in a sample of Class I sources. 

5.1 .4. Correlations with jet/outflow parameters 

Forbidden atomic lines of [O I] (6300 A) and [Fe II] (1.64 /urn) 
are used to trace t he outflow activity in youn g stars (e.g., 
Cabr it et alJ 1 19901: lAgra-Amboage et all 120091) . iGiidel et alJ 
(1201 Oh investigated this correlation to test the rela t ion be tween 
[Ne II] emission and outflow activity. Shan g et all (HpiO) mod- 
eled the flux and shape of the [Ne II] and [O I] lines and found a 
clear correlation between both lines, indicating that neon would 
also trace jets. 

We have tested the correlation between the luminosity of the 
[Ne II] line and the luminosity of the [O I] line at 6300 A com- 
piled from the literature, obtaining a positive result. We have also 
te sted a correlation w ith the product of Lx and L[o i] , as proposed 
in Gud el etalJ (l2010). but a correlation was not found. We have 



further used the luminosities of the [Fe II] line at 25.99 /-im ob- 
tained in our sample to study the dependence of [Ne II] on this 
tracer, and no correlation was found. We emphasize, however, 
that the tests were done with a large number of upper limits due 
to the lack of detections in the tested parameters. 

Three objects with detected [Ne II] emission also show 
[Fe II]: L1551 IRS 5 (Class I), DG Tau (Class I/II), and 
SST 041412.2+280837. Since [Fe II] is a known tracer of 
shocks, the origin of the [Ne II] emission in these particular cases 
might be explained by the shock scenario. However, this is not 
clear for the rest of the objects with [Ne II] emission for which 
there is no detected emission from [Fe II]. 



5.1.5. Correlations with L x xM wind 

According to the model presented in lShang et alJ (l2010h . [Ne II] 
is produced in jets from young stars. Based on the X-wind model 
theory, they obtained line profiles for [Ne II], [Ne III], and [O I]. 
They tested the dependency of the line luminosities on critical 
parameters of this model such as Lx and the wind mass loss 
rate (M w i n d)- The best correlation was found for the product of 
both parameters LxX M w j n( j. The relation between wind mass 
loss rate and accre tion rate is M w ; n d = fx M acc , with / w 0.25 
(IShang et al.ll2004h . However, in our sample, no such correlation 
is found. 
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5.2. [Fe II] emission 

[Fe II] is a tracer of shock activity typically observed 
in the near-IR at 1.53 and 1.64 fan (iDavis et al.1 l2003t 
lHartigan & Hillenbrand! 12009b . Its excitation mechanism can 
also be stimulated by the irradiation of protoplan etary disks by 
high-energy photons dGorti & Hollenbachll2008l) . We have de- 
tected [Fe II] emission in 7 sources. In five cases, we have de- 
tected both iron lines (at 17.93 and 25.99 /mi). Typically, when- 
ever both lines have been detected, the luminosity of the [Fe II] 
line at 25.99 /mi is roughly a factor of 3 higher than the lumi- 
nosity of the line at 1 7.93 A mi. This result is consistent with 
IHollenbach & McKed (fl989). where in a shock environment the 
[Fe II] line at 25.99 /mi is a more effective coolant than the 
line at 17.93 /mi. Four out of the seven sources with detected 
[Fe II] lines (L1551 IRS5, DG Tau, RW Aur, IRAS 04239+2436, 
and IRAS 04303+2240) have previous indications of jet activity 
based on near-IR and/or optical observations (see Appendix lAl. 
We have tested the correlation between the luminosity of the 
[Fe II] line (25.99 /mi) obtained in this study with different pa- 
rameters: stellar mass, luminosity, X-ray luminosity, infrared lu- 
minosity, mass accretion rate, and luminosity of the [O I] line 
(see Table 0. In some cases (mass and luminosity) the results 
from the different correlation tests for L[p e ii] give largely dif- 
ferent results. We remind that mostly upper limits are available 
for the [Fe II] line; indeed, we have only 7 detections for 53 
to 64 points, depending on the parameter tested. Furthermore, 
most correlations were found with the Spearman method, but 
the Kendall and Cox methods are in general more robust and no 
correlations were found with these methods. In summary, except 
for mass accretion rate, no correlation is found for the [Fe II] 
line luminosity. 



5.3. H 2 emission 

Three pure rotational lines of Fb are covered by the Spitzer IRS: 
12.28 /mi (S2), 17.03 /an (SI), and 28.22 /mi (SO). Detections 
of H2 in the near-IR (S(l) tra nsition at 2.12 / i rn) were reported 
in a handful of objects (e.g. Iltoh et al.ll2003t lBar"v et al.ll2002i 
I2008bt ICarmona et al]|2008at1Panic et al.ll2009l) . In addition to 
H2 disk emission, H2 emission may also originate in shock envi- 
ronments; in particular the near-IR lines and the mid-IR line at 
12.28 nm are associated with shock densities of n~ 10 4 ~ 5 cm -3 
dHollenbach et al.ll!989b IHollenbach & McKedll989h . 

We have detected H2 in six sources: one of the newly 
identified Taurus members SST 043905.2+233745, which is a 
Class I source, the Class II sources FS Tau A, FX Tau, IQ Tau, 
V710 Tau, and the Class III HBC 366. FX Tau and V710 Tau are 
both binary systems composed of a classical and a weak-lined 
T Tauri star. These are the first detections of H2 reported toward 
these sources. 

We have tested correlations between the total luminosity of 
the H2 line with the following parameters: L*, M*, Lir, Lx, M aC c, 
and L[o i] . The total luminosity of the H2 line was calculated by 
adding the luminosities of the three lines accessible in the IRS 
wavelength range. Our analysis does not reveal any correlation 
with the above parameters. Although correlations with H2 were 
found for both the mid-IR and [O I] luminosities, this was only 
the case for one test method out of three. We prefer, thus, to 
conclude that no correlation was found. 



6. Summary and Conclusions 

We have searched for gas in a large sample of pre-main sequence 
stars (64 objects) spanning a broad range of evolutionary stages. 
Our sample is dominated by Class II objects but it also includes 
Class I, III, and Herbig AeBe stars. We have detected 36 gas 
emission lines; in particular, we have obtained 18 detections of 
the [Ne II] line. Nine out of the 18 [Ne II] detections are reported 
here for the first time. For 4 targets there are previously reported 
detections of [Fe II] based on near-IR observations, which we 
complement with the mid-IR information. The six sources in 
which we detect H2 emission do not have previously reported 
detections. 

The [Ne II] detections are concentrated towards Class II 
sources; we have not detected [Ne II] lines towards Class III or 
Herbig AeBe stars. Two Class I sources show [Ne II] line emis- 
sion. These sources have higher accretion rates, and are more 
likely to show jet/outflow phenomena, but they also are more 
embedded in the circumstellar material. If neon line emission is 
produced by jets, we would expect to have a larger detection rate 
for Class I objects. We have attempted to disentangle sources 
with and without jets, but we cannot statistically separate the 
two populations. 

We have tested correlations between the [Ne II] line with 
parameters probing different stellar-disk properties; high-energy 
photon irradiation, accretion, and jets/outflows. We have found 
that the luminosity of the [Ne II] line is dependent on the con- 
tinuum mid-IR lum inosity. The dependence with Lx obtained in 
iGtidel et all (1201 Ol) is confirmed for the sub-sample of Class II 
sources. A go od correlation is als o found for the L[o rj, echoing 
the finding by Gud el et aD d2010l) and clearly supporting a trend 
with jet parameters. 

In the jet/outflow scenario proposed by 
IHollenbach & McKeel d!989l) . the [Ne II] line is proposed as 
an effective coolant. For high density shocks (n~ 10 5 " 6 cm" 3 ), 
the [Fe II] line (25.99 //m) is more effective, with expected 
luminosities up to 2 orders of magnitude higher. Only in lower 
density environments (n~ 10 3 " 4 cm" 3 ) and high velocities 
(> 80 km s" 1 ) the luminosity of the [Ne II] surpasses the lumi- 
nosity of [Fe II]. Among the sources with [Ne II] detections, 
three also show [Fe II] emission: DG Tau, L1551 IRS 5 (both 
known to drive jets), and SST 041412+289837 (for which 
little information is available in the literature). In these cases, a 
plausible explanation for the [Ne II] emission is outflow activity. 
The question arises for sources with [Ne II] detections that 
do not show [Fe II] emission: if [Ne II] is excited in shocks 
produced by outflows, what is the explanation for the absence of 
[Fe II] emission? This might be due to different reasons; perhaps 
lower signal-to-noise was obtained in the [Fe II] wavelength 
range, or the conditions in the emitting region are such that the 
neon line is brighter, or simply for those sources, neon is not 
excited by shocks, but its line emission originates, e.g., from the 
disk atmosphere. 

The [Fe II] emission line has been detected in a smaller num- 
ber of sources, in comparison to the [Ne II] line. The statistical 
analysis has shown that the luminosity of [Fe II] is dependent on 
the mass accretion rate. This indicates a link with outflow phe- 
nomena, since the outflow rate is proportional to the mass accre- 
tion rate. On the other hand, no correlation between the [Fe II] 
luminosity and the luminosity of the outflow-tracing [O I] (6300 
A) line was found. 

Furthermore, Spitzer IRS can access three H2 emission lines. 
We have obtained few detections, and, despite the large number 
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Table 2. Correlations tested for the luminosity of the [Fe II] line. 



All infrared classes 

N° Upper L. Probability 



Var. 


N° Points 


Total 


[Fell] 


Cox 


Kendall 


Spearman 


M, 


53 


48 


48 


0.55 (0.45) 


0.17(1.3) 


0.001 (0.44) 


L, 


57 


51 


51 


0.77 (0.08) 


0.14(1.5) 


0.001 (0.46) 


Lx 


47 


43 


40 




0.66 (0.44) 


0.25 (0.17) 


Lir 


64 


57 


57 


0.81 (0.05) 


0.62 (0.49) 


0.49 (0.09) 


M acc 


47 


42 


39 




0.01 (2.6) 


0.0004 (0.52) 


L[oi] 


27 


24 


24 




0.16(1.4) 


0.09 (0.33) 



Notes. Probabilities (P) of a correlation achieved by chance are presented in columns 5, 6, and 7, where P < 5% are shown in bold fonts. The value 
of the correlation parameters are shown between parentheses. 

Table 3. Correlations tested for the luminosity of the H2 line. 



All infrared classes 







N° Upper L. 




Probability 




Var. 


N° Points 


Total 


H 2 


Cox 


Kendall 


Spearman 


M. 


53 


48 


48 


0.97 (0.001) 


0.79 (-0.26) 


0.57 (-0.078) 


u 


57 


52 


52 


0.17 (1.9) 


0.61 (-0.51) 


0.41 (-0.11) 


Lx 


47 


42 


39 




0.86 (-0.17) 


0.54 (-0.095) 


Lir 


64 


58 


58 


0.03 (4.5) 


0.09 (1.7) 


0.09 (0.01) 


Macc 


47 


44 


40 




0.22 (-1.2) 


0.43 (-0.15) 


L[oi] 


27 


26 


26 




0.99 (0.00) 


0.04 (0.40) 



Notes. Probabilities (P) of a correlation achieved by chance are presented in columns 5, 6, and 7, where P < 5% are shown in bold fonts. The value 
of the correlation parameters are shown between parentheses. 



of upper limits, we have tested different correlations. Our results 
do not show any dependency on the parameters tested. 

Our results show a trend towards accretion-related phenom- 
ena as the origin for the forbidden lines. Several mechanisms 
could be playing a role in the heating and excitation mechanism 
of gas in protoplanetary disks. A sample with several Class I, 
and Class III objects would be more suitable to study possible 
trends with evolutionary state. But no matter what mechanism 
produces the neon emission, it is likely originating very close 
to the star. High spatial and spectral resolution observations are 
crucial in order to determine its origin. 
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Fig. 3. Sources with [Fe II] detections at 17.94 and 25.99 /urn. We recall that the spectrum of SST 041412.2+280837 was obtained 
using the low resolution modules (SL and LL). 



12 



Baldovin-Saavedra et al.: Searching for gas lines in the young population of the TMC 



io J1 fi as 1 1 

* Class II, 

• Class I/II 
T Class III 
■ HAe 



10" 



10 a 



10" 




10 J 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 
M, [M slm ] 



10" 



10 B 



10" 



A 'Class' 1 ' 
* Class II 



10 J 



10- 5 ' 




10" 10-* 
L, [erg/s] 



A Class I 

★ Class II 

• Class I/II 
T Class III 
■ HAe 




io- 5 



10 J1 



10" 



10"' 



: A Class '1" 
: ★ Class II 



10" 



10"' 




10"" 10 30 1 31 
L x [erg/s] 



10 J 



10 s ' 



10"' 



10" 



AClass I 

★Class II 

• Class I/II 

TClass III 
■ HAe 



i V 



10" 



10~ 9 10~ 8 10~ 7 10" 

Mace [ M sun/y r ] 



10" 



io- 5 



10 J1 



10" 



10" 



A Class I 

★ Class II 

• Class I/II 
T Class III 
■ HAe 



10" 



0.1 1.0 10.0 100.0 

Ha EW [A] 



Iff" E A 'Class '("' 
; ★ Class II 
• Class I/II 
T Class III 
. ■ HAe 

io 30 . 



10" 



10" 




io" 7 10 ZB 10" 9 10 
L[Fe II][erg/s] 



30 . n 31 



10 J 



10 al 



10" 



10" 



Class I 
Class II 
Class I/II 
Class III 
HAe 



it i 



Mo i] 



[erg/s] 2 



10 J 



10 



: AClass I 

: ★Class II 

' •Class I/II 

▼ Class III 
■ ■ HAe 



* 
* 



10" 



10" 



* ** 

lif 



10" u 10"' 10" 10" 10" 10" 
L x * M.i„d [erg/s M aun /yr] 
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Fig. 5. Same as Fig.|4]for Class II (stars) and sources classified as I and II in the literature (circles). 
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Table 1. Stellar properties used in this study. 
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2.41 


8.90 


41 


-7.24 


8.05 


59.5 


-2.63 


V 


1,22,27 


UX Tau A 




1.10 


1.81 


9.5 


-9.30 


0.95 


1.27 


< -4.97 


X 


2,6,9,22 


UZTau 




0.47 


0.49 


73 


-8.01 


0.89 


3.98 




X 


1 


VY Tau 


III 


0.57 


0.55 


4.9 




1.29 


0.38 


< -5.56 


X 


2,7,10,22 


XZ Tau (A,B) 




0.34 


0.33 


62 


-8.74 


0.96 


13.1 




V 


1,2,25 


CokuTau-1 




0.70 


0.15 


70 


-7.36 




3.90 




V 


1,31 


Coku Tau-3 




0.46 


0.98 


5 




5.85 


1.31 




X 


1 


Coku Tau-4 




0.43 


0.68 


1.8 




< 0.21 


0.70 




X 


2,9,10 


AB Aur 


HAe 2.40 


49.00 33 


-6.85 


0.35 


64.6 


-3.36 


V 


1,11,28 


GM Aur 




0.52 


0.74 


96 


-8.02 


1.68 


0.68 


-4.90 


X 


3,5,7,22 


RW Aur A 




1.34 


3.37 


76 


-7.50 


1.60 


4.70 


-3.11 


V 


4,23,25,32 


SU Aur 




1.91 


9.90 


63 


-8.20 


9.46 


15.0 


< -3.85 


X 


1,22 


UY Aur (A,B) 




0.57 


1.25 


36 


-7.37 


< 0.25 


10.6 


-3.62 


V 


2,9,12,23,25 


Haro 6-28 




0.21 


0.12 


48-92 


-8.10 


0.25 


0.47 




X 


1,12 


Haro 6-37 




0.91 


1.74 


22 


-8.12 




2.35 


< -5.07 


X 


2,6,22 


L1551 IRS5 




1.58 


2.60 


83 


-6.47 


0.021 


62.0 




V 


1,13,14,30 


MHO- 1/2 


I/II 


0.34 


0.42 


58 


< -8.48 


1.64 


7.20 




X 


1 


V710Tau 




0.68 


1.10 


69 




1.38 


0.75 




X 


1,6 


V773 Tau 




0.58 


5.60 


2 


< -8.6 


9.49 


4.46 


-3.93 


X 


1,6,22 


V892 Tau 


HAe 2.89 


77.00 7-13 




9.21 


132 




V 


1,29 


HBC 366 


III 


0.40 


2.60 


1.3 


< -8.94 


4.14 


0.23 




X 


1,16 


HD 31648 


HAe 2.30 


23.7 


18.3 


-8.39 




34.8 




X 


17 


HD 35187 


HAe 2.30 


1.37 


5.1 


-7.81 




18.5 




V 


11,18 


HD 36112 


HAe 2.20 


1.45 


14.5 


-6.70 




17.4 




X 


18,19 
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Source 


Class Mass 
(M ) 


L 


EW(Hq-) log(M/M yr~') 
(A) 


Ly 

(iu erg/s; 


1 m lA TD 

^-rnia— IK 

no 32 

erg/s) 


lo2(Lm n/Ln) Jet a 


Ref. 


IC2087IR 


II 






22 


-6.20 


4.67 


16.3 


X 


1 


IRAS 04016+2610 


I 


1.45 


3.70 


41-56 


-7.15 


4.47 


13.6 


V 


1,21,27 


IRAS 04101+3103 


HAe 






11 






10.7 


X 


15 


IRAS 04108+2803 B 


I 


0.36 


0.4 






0.42 


2.46 


V 


1,27 


IRAS 04239+2436 


I 




1.27 








8.76 


V 


13,27 


IRAS 04278+2253 


I 




7.20 


17 


-5.18 




28.6 


X 


4,20,21 


IRAS 04299+2915 


II 






6.5 






0.47 


X 


15 


IRAS 04303+2240 


II 


0.51 


2.20 


67 


-6.05 


5.01 


6.43 


X 


1 


IRAS 04361+2547 


I 




3.70 






<0.30 


8.28 


v 


1,9,27 


IRAS 04365+2535 


I 




2.20 






1.99 


6.13 


v 


1,27 


SST 041412.2+280837 


I 












1.73 


X 


33 


SST 042936.0+243555 


II 












0.32 


X 


24,33 


SST 043905.2+233745 


I/II 












0.23 


X 


24,33 



V indicates presence of jet, x indicates that no jet has been reported. 
(i) indicates no information available. 



Refer ences. rilGudel et alJ d2007h : [2] IWatson et al.1 ( 120091) : [3] iGullbring et all dl998h: f41 IWhite & Hillenbrand! d2004h: T51 lHerbig& Belli 
d!995h: T61 IWhite & Ghe3 d2001h: [7] Gudel priv ate communicati on 2009; [81 [Nguyen et alj d2009h: [91 iNeuhauser et al.1 d!995h; [10 ] 
Cohe n & KuhiHl979h: ri ll barcia Lopez et aljd2006h: [12 1 compiled in M cCabe et al.1 d2006h ; [131lPrato et all f"2 009h 114l lKenvon et al.l dl998h: 
T151 iKenvon et al.l d!990l); [16] comp iled in iMartin et al.1 d!994h: fl71 Mannings & Sargend fT997l); [181 iManoi et al.1 d2"006h: fl9 1 llsella et al.1 
d2008h: T201 iBitner et al] d2008h: T211 iFurlan et al] d2008h. T221 ICabrit et al] d!990). 1231 Haiti gaii et al] dl995h, [ 241 iLuhman et al] d2010l); [25 1 
iHirth et al] dl997h: T26 1 lHartigan et al] d2004h: [ 271 iMoriarty-Schieven et al] dl9 92). [28] Corcoran & Ravi dl997h . [291 iMaheswar et al] d2002). 
[30] Davis etalj OOOl) . [31] Eisloffel & Mundtl dl998h . [321 lGiidel et al] d2010h . r331 lRebull et al] d2O10h 
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Table 2. Line luminosities derived in this study. 





H 2 


[Ne II] 


[Ne III] 


H 2 


[Fe II] 


[S III] 


[Fe II] 


H 2 


[S III] 


Name 


12.2786 


12.8135 


15.5551 


17.0348 


17.9359 


18.7130 


25.9883 


28.2188 


33.4810 


AATau 


1.2 


2.7 ± 0.3 


2.2 


2.8 


1.7 


2.2 


1.7 


1.9 


1.9 


BPTau 


2.1 


1.0 


1.3 


3.5 


2.3 


1.8 


4.8 ± 0.4 


1.9 


1.9 


CWTau 


11.2 


4.9 


9.5 


8.2 


10.9 


4.9 


3.2 


6.2 


8.8 


CZ Tau (A,B) 


4.2 


3.8 


3.0 


3.9 


6.5 


5.4 


2.8 


3.2 


3.6 


DF Tau (A,B) 


5.2 


8.4 


2.4 


10.0 


3.6 


4.6 


2.3 


3.0 


5.8 


DG Tau 


13.6 


43.0 ± 4 


9.5 


9.8 


10.3 


10.6 


40 ±2 


15.2 


9.7 


DG Tau B 


5.1 


10 ±2 


12.3 


11.6 


8.5 


10.1 


6.1 


3.8 


12.9 


DK Tau 


5.7 


5.6 


6.7 


7.9 


6.5 


2.7 


3.3 


3.1 


7.4 


DL Tau 


7.7 


4.0 


5.5 


5.9 


6.6 


5.6 


4.7 


2.0 


8.9 


DM Tau 


1.1 


0.8 ± 0.2 


1.2 


0.8 


1.2 


1.0 


0.6 


1.5 


0.3 


DN Tau 


1.3 


1.2 


0.8 


0.9 


0.8 


0.9 


1.9 


1.4 


2.0 


DO Tau 


4.6 


4.5 


3.7 


4.3 


3.2 


4.6 


2.2 


3.9 


8.0 


DR Tau 


12.6 


6.5 


12.9 


18.0 


8.8 


8.8 


3.7 


3.3 


14.2 


FS Tau A 


10 ±2 


17 + 2 


4.9 


5.2 


7.9 


7.0 


2.6 


4.1 


4.0 


FV Tau (A,B) 


6.7 


4.3 


4.2 


6.6 


6.3 


1.2 


3.6 


2.1 


5.4 


FX Tau 


1.5 


1.1 


3.8 


2.5 


3.3 


5.3 


1.8 


2.3 ± 0.4 


1.2 


FZ Tau 


8.2 


4.7 


5.2 


4.7 


8.0 


7.0 


2.6 


6.6 


5.8 


GG Tau 


7.0 


6.5 


4.1 


8.2 


5.1 


3.9 


3.2 


3.8 


8.9 


GI Tau 


2.0 


2.0 


2.9 


3.2 


2.8 


2.1 


0.9 


2.2 


3.1 


GK Tau 


3.0 


5.0 


2.5 


2.8 


3.4 


2.1 


2.3 


3.9 


2.9 


GV Tau (A,B) 


50.5 


40.0 


41.6 


21.4 


76.4 


30.2 


39.7 


56.8 


51.6 


HL Tau 


40.0 


14.1 


14.5 


19.8 


14.7 


13.9 


22.8 


17.1 


35.2 


HN Tau 


6.3 


2.1 


4.1 


3.3 


5.6 


4.0 


3.4 


2.3 


6.0 


HP Tau 


6.8 


4.6 


4.8 


3.9 


8.8 


6.5 


5.0 


2.6 


6.9 


HQ Tau 


5.4 


5.2 


6.4 


5.9 


10.7 


3.1 


4.0 


3.0 


8.4 


IQ Tau 


1.5 


1.2 ± 0.4 


1.7 


1.6 


1.4 


2.5 


1.6 


3.9 ± 0.6 


1.7 


RY Tau 


36.6 


13.7 


33.7 


25.2 


17.0 


23.3 


7.9 


7.7 


40.9 


T Tau(N,S) 


49.5 


400 ±9 


21.9 


54.6 


31.2 


25.4 


36.9 


109 


119 


UX Tau A 


6.9 


5.6 


2.5 


2.9 


7.5 


7.4 


1.8 


1.8 


1.3 


UZTau 


9.1 


5.1 


4.5 


7.3 


8.6 


9.1 


4.5 


9.7 


8.8 


VY Tau 


0.3 


0.4 


0.3 


0.2 


0.5 


0.5 


1.7 


1.3 


1.7 


XZ Tau (A,B) 


9.5 


7.2 ± 0.9 


5.4 


9.2 


15.9 


9.9 


4.6 


9.6 


9.6 


Coku Tau-1 


4.1 


27+1 


3.9 


3.7 


5.7 


5.2 


5.1 


7.1 


6.0 


Coku Tau-3 


1.5 


1.5 


0.5 


0.8 


1.5 


1.3 


1.9 


1.6 


1.3 


Coku Tau-4 


1.3 


3.2 ± 0.2 


0.8 


0.9 


0.9 


1.6 


1.4 


1.1 


2.1 


AB Aur 


29.7 


16.8 


18.5 


25.1 


18.2 


34.9 


44.3 


25.7 


51.8 


GM Aur 


2.6 


1.5 ± 0.3 


1.0 


1.0 


1.4 


2.2 


1.3 


1.1 


4.0 


RW Aur 


14.6 


3.8 


7.9 


12.9 


10 ±2 


6.6 


29 ± 1 


8.1 


5.0 


SU Aur 


25.2 


6.5 


6.8 


7.8 


18.6 


17.4 


9.8 


11.2 


10.3 


UY Aur (A,B) 


14.1 


11 ±2 


9.9 


8.7 


10.0 


5.4 


4.3 


8.3 


13 


Haro 6-28 


4.4 


4.4 


10.5 


3.7 


3.7 


5.5 


1.6 


1.4 


3.9 


Haro 6-37 


5.4 


4.0 


2.2 


4.9 


5.5 


5.7 


2.1 


3.0 


4.0 


L1551 IRS5 


30.4 


129 ±9 


75.8 


23.3 


200 ± 19 


49.8 


597 ±22 


20.1 


13.1 


MHO- 1/2 


8.4 


10 ±2 


4.7 


6.9 


19.7 


12.8 


3.5 


7.7 


11.3 


V710 Tau A 


0.9 


0.4 


0.6 


1.0 ± 0.3 


0.5 


1.1 


1.0 


2.1± 0.4 


3.5 


V773 Tau 


3.6 


34 ±2 


4.7 


8.9 


6.8 


4.1 


3.3 


3.6 


6.4 


V892 Tau 


100.4 


40.1 


141.4 


53.7 


79.9 


56.9 


84.3 


64.3 


166.8 


HBC 366 


0.4 


0.2 


0.1 


0.5 ± 0.1 


0.3 


0.2 


0.3 


0.5 


1.0 


HD 31648 


45.3 


7.5 


9.9 


15.4 


35.1 


15.5 


7.1 


5.3 


10.4 


HD 35187 


11.6 


11.4 


10.8 


9.7 


16.9 


13.5 


8.4 


8.3 


11.1 


HD 36112 


21.0 


11.1 


11.2 


7.7 


15.8 


18.3 


9.5 


5.4 


1.8 


IC2087IR 


20.9 


11.3 


16.3 


11.8 


7.8 


7.2 


4.5 


6.4 


5.9 


IRAS 04016+2610 


9.7 


5.6 


15.9 


9.1 


8.5 


11.4 


15.3 


8.0 


21.1 


IRAS 04101+3103 


4.8 


6.3 


6.1 


3.8 


9.5 


7.1 


6.2 


8.7 


4.9 


IRAS 04108+2803 B 


3.1 


3.9 


4.6 


4.6 


6.5 


4.5 


8.0 


6.6 


10.0 


IRAS 04239+2436 


4.6 


4.4 


6.0 


7.3 


24 ±3 


8.3 


83 ±3 


16.4 


38.8 




1 A 1 

L<\.L 


1-7 1 
1 /. i 


1 O 1 

ty.t 


C -7 
J. / 


Tin 
21. y 


1 1 A 

13.4 


D.U 


n -7 

y. 1 


1 -7 O 

1 1.0 


IRAS 04299+2915 


3.6 


3.2 


3.3 


3.6 


4.6 


2.0 


2.9 


2.7 


6.4 


IRAS 04303+2240 


14.8 


8.2 


11.5 


13.9 


21 + 2 


9.8 


48 ± 1 


16.9 


11.1 


IRAS 04361+2547 


7.5 


3.0 


7.8 


7.0 


8.9 


6.6 


6.8 


34.3 


26.4 


IRAS 04365+2535 


9.4 


4.6 


8.4 


3.7 


10.5 


6.5 


11.9 


12.5 


19.0 


SST 041412.2+280837 


4.0 


7.5 ± 0.6 


7.1 


3.2 


5.3 ± 0.7 


2.8 


3.5 ± 0.5 


4.7 


6.6 


SST 042936.0+243555 


0.9 


2.5 ± 0.2 


0.7 


0.8 


0.5 


0.8 


0.6 


0.8 


1.0 


SST 043905.2+233745 


0.8±0.2 


1.6 ± 0.1 


0.4 


1.2 ± 0.1 


0.8 


0.09 


0.5 


1.9 


2.6 



Notes. Detections plus lcr errors are in bold characters, 3tr upper limits are presented in normal font. All values are given in units of 10 28 erg s 
The wavelength of the emission lines are given in /im. 
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Appendix A: Results for Individual Objects 

We start with L1551 IRS 5 because it is a good example to illus- 
trate the procedure used to study the line emission. The rest of 
the objects are ordered as in Table Q] 

L1551 IRS 5 is a Class I protostar driving a molecular outflow. 
High angular resolution observ ations in the radio have re vealed a 
binary with separation of 0".3 dBieging & Cohenlll985l) . One of 
the binary compone nts is now known to be a binary itself with 
separation of 0".09 (Lim & Takakuwa 2006). We do not have 
any further information on the nature of the components of the 
system; Spitzer observations do not resolve the components (its 
spatial resolution is ~ 4" at 13 fim). 

This target is part of a GTO program. The observations were 
done with 6 different pointings offset in spatial direction. The 
pointings correspond to three sets of standard nod observations. 
Therefore we could obtain three average spectra located at three 
positions; one centered at the star position, and the two oth- 
ers at the offset positions. We can use these data to compare 
the continuum-subtracted flux at the different positions in order 
to have a rough estimate for the contribution of the extended 
emission, whenever the line is detected at multiple positions. In 
Fig. lA.ll we show the spectrum of L1551 IRS 5 around the po- 
sition of the [Ne II] line. Each plot in the left panel represents 
the average of the spectra between 2 nods at a given position. 
In this specific case, we have three different positions identified 
as "Offset 1", "Offset 2", and "Centered". We have labeled each 
position in the upper left side of each plot. The maximum flux in 
the continuum is achieved for the "Centered" position. The line 
is detected in all three sets of pointings and we point out that it is 
shifted by 0. 1 pm with respect to the expected wavelength of the 
[Ne II] line. Such a wavelength shift corresponds to a velocity 
shift of ~ 230 km s -1 . This result is co nsistent with the ve locity 
of the [Fe II] jet (285 km s~' ) de r ived b y|DavisetaD d2003l) from 
near-IR observations. iPvo et al.l(l2009l) have observed the [Fe II] 
jet in the near-IR and derived a position angle (P. A.) of 260° 
for the northern component and 235° for the southern compo- 
nent. Spitzer observations have been obtained using a slit R A. 
of -56.7° and -141.5° for the SH and LH modules respectively. 

In the right panel of Fig. lATl we have plotted the continuum- 
subtracted flux. In this way we could verify that although there 
is an extended component in the emission, the maximum, i. e., 
the highest line peak, corresponds to the position centered on the 
object. But we recall that extended emission close to the star can 
be present and not be resolved by Spitzer. 

In the case of [Fe II] at 17.9 and 25.9 pm, the lines are de- 
tected in all three pointings (Fig. IA.21 . Both lines are centered 
at the expected wa velength, in contrast to near-IR observations 
dDavis et al.ll2003l) and the [Ne II] line. In both cases, there is 
an important contribution from extended emission inferred from 
the high peak of the continuum-subtracted line in the position 
labeled Offset 2, suggesting diffuse emission as the main con- 
tributor to the line luminosity. However, the highest luminosity 
is still obtained from the pointing centered at the position of the 
star 

iDavis et al ] (120031) have reported a detection of an H2 emis- 
sion line at 2. 12 /mi, but we do not detect H2 in the mid-IR spec- 
tra obtained with Spitzer. 

The emission lines from L1551 IRS 5 are the brightest of the 
sample, with luminosities on the order of 10 30 erg s _1 . We have 
obtained VISIR/VLT observations in order to study the kinemat- 
ics and profile of the [Ne II] line. We have not detected [Ne II] 
emission at the position of the star, which might suggest a jet ori- 
gin for the [Ne II] emission, although we have detected the sys- 
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Fig. A.l. Spectrum of L1551 IRS 5 surrounding the position of 
the [Ne II] line at 12.8 pm. The plots on the left panel show the 
average of 2 pointings at each position. We have labeled the posi- 
tions as Offset 1, Centered, and Offset 2. The spectrum obtained 
for the position "Centered" has the highest flux in the continuum 
so it was used for the line analysis. The expected position of the 
[Ne II] line is overplotted with a black dash-dotted line, we note 
that the line is shifted. The right panel shows the difference be- 
tween the line flux and the continuum flux in order to show the 
real strength of the line. 



tem in imaging mode using the [Ne II] filter (Baldovin-Saavedra 
et al. 201 1, in preparation). 

AA Tau is a Class II source for which [Ne II] i s the only line 
we de tected (of the lines we studied). However, ICarr & Najital 
(2008) have presented the spectrum of AA Tau showing multi- 
ple emission lines from organic molecular species (HCN, C2H2, 
and CO2), wat er vapor, OH, and [Ne II]. Follow-up observations 
(at R~ 80 OOO. lNaiita et al.l2009l) have confirmed that the [Ne II] 
emission is centered at the stellar velocity, which means that its 
origin lies in the disk rather than in a jet or outflow. Nevertheless, 
due to the high inclination of the system (~ 75°, iBouvier et al.l 
1999), the absence of velocity shifts i n the emission lin e does not 
rule out a jet origin for the emission. ICox et al" I (12001 have in- 
deed reported a microjet detected using coronographic imaging 
with the Hubble Space Telescope (HST). 

BP Tau is a Classical T Tauri star (CTTS) that was observed 
together with background observations. We detect [Fe II] at 
25.99 pm, with a marginal detection at 17.93 pm. There is no 
record in the literature of previous [Fe II] detections. We have 
plotted in Fig. IA.3I the SH spectra of AA Tau together with 
BP Tau and DK Tau. Both spectra show spectral f eatures similar 
to the ones observed in AA Tau and reported by ICarr & Najital 
(120081) . The presence of molecular lin e emission in the spectra o f 
young stars were recently reported in lPontoppidan et al] (|2010). 
Although the spectra of BP Tau and DK Tau are quite similar to 
the spectrum of AA Tau, we do not detect [Ne II] emission from 
these two stars. iGtidel et al.1 (1201 Oh detected [Ne II] emission in 
BP Tau using the same data set as this study, although the line is 
weak (~ 10 27 erg s~'). We observe an excess at the position of 
the [Ne II] line (see Fig. [A3}, but below our 3cr detection thresh- 
old. A similar situation occurs for the H2 line at 12.27 pm (see 
Fig. [A3}. The presence of numerous faint emission lines likely 
makes our detection threshold slightly too high. 
DG Tau is classified as both Class I and II in the literature. This 
young star drives a collimated Herbig-Haro (HH) jet observed 
in the near-IR. The Spitzer observations were made in the same 
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Fig.A.2. Spectrum of LI 551 IRS 5 surrounding the position of 
the [Fe II] line at 17.9 and 25.9 pm. The right panel shows the 
difference between the line flux and the continuum flux in order 
to show the real strength of the line. The plots on the left panel 
show the average of 2 pointings at each position. Labels follow 
the convention used for Fig. IA.ll 



way as L1551 IRS5, with 3 sets of 2 nod pointings. The [Ne II] 
line is detected in all of them. The observations were made us- 
ing a slit P.A. of -57.1° and -141.9° for the SH and LH mod - 
ules respectively. The jet P. A. is -138° dLavallev et al.lll997b . 
almost coincident with the LH module P. A. When comparing 
the difference between the line and the continuum flux for each 
set of pointings, we note that the bulk of the emission comes 
from the position centered on the coordinates of the star, with 
the ratio being roughly 2:1 between the centered and the offset 
positions. A detect ion of [Ne II] has been previously reported in 
Giide l et al] (|2010), based on a different data set than this study. 
We detect [Fe II] emission at 25.99 pm; this line is detected in 
all pointings, again with a higher contribution from the position 
centered on the coordinates of the star. Emissi on from [Fe II] 
in the near-IR has been previously rep orted (e.g.. | Bacciotti et all 
120021: iDavis et al.ll2003l) . In particular, iDavis et al.l d2003l) do not 
detect any [Fe II] emission within 1" from the source, suggesting 
a jet origin for the [Fe II] emission. 

DG Tau B has an edge-on disk and a bipolar jet detect ed in 
the op tical dMundt & Friedll 19831: lEisloffel & Mundtll 19981) and 
radio dRodripuez et al.l fl995). The jet is oriented perpendicular 
to the disk. The IRS spectrum has been obtained at only one 
position. We have detected [Ne II] emission towards this object. 
DM Tau is a Class II object with [Ne II] line emission. The 
IRS observations were done with a background observation, so 
we are confident that the detected line comes from the source it- 



self. Based on Spitzer observations made wit h the low-resolution 
modules (SL and LL), Espai llat et al.l (120071) have derived a lu- 
minosity of the [Ne II] line of 1.3 x 10 28 erg s _1 . In contrast, 
we have obtained a luminosity of (8.0 ± 2.0) x 10 27 erg s~', 
which is 40% lower than what is presented in Espailla t et al.l 
(2007). This object is particularly inter esting because it belongs 
to the category of transitional disks (ICalvet et al.ll2002l) . These 
objects show a lack of emission at wavelengths shorter than 10 
pm, which is sometimes interpreted as an evidence of an inner 

hole and planet growth. 

FS Ta u A is a close binary with separation of 0".25 (Rris t et al.l 
1998). The system is surrounded by a filamentary nebula. Spitzer 
observations have been done in 3 sets of 2 nod pointing; the 
[Ne II] line is detected in all of them. By analyzing the flux 
difference between the continuum and the line, we have found 
that the height of the line in the centered position is comparable 
to the height of the line at one of the offset positions, indicat- 
ing the emission is likely extended. In addition we have detected 
H2 at 12.28 pm, while there is no detection of the lines at 17.03 
and 28.22 pm. At 28.22 pm, there could be a marginal detec- 
tion, but only below the detection threshold of 3<x. In the case of 
the 12.28 pm line, it has been detected in 2 out of the 3 sets of 
pointings. Comparing the continuum-subtracted line flux at the 
different positions, we determine that the line flux is higher for 
a position shifted with respect to the star position, indicating an 
important contribution from extended emission. 
FX Tau is a binary with 0."89 separati on composed of a 
classical and a weak-lined T Tauri star (IWhite & Ghezll200l . 
Assuming that the mid-IR emission is dominated by the CTTS, 
we have included this system in the Class II objects for our anal- 
ysis. We have detected H2 emission at 28 pm towards this object, 
and it is the only line detected. No previous detections have been 
reported. 

IQ Tau is a CTTS showing [Ne II] line emission. There is no 
information in the literature on possible jet or outflow emission 
associated with this object. Carmona et al. 201 1 (in preparation) 
detected [O I] emission towards this object employing optical 
spectroscopy (R~ 40 000), observing a velocity shift of the line 
which indicates the presence of an outflow. 
T Tau is a triple system visible in the optical and bright in 
X-rays. T Tau S is an infrared companion located 0".6 away 
(lDvcketal.1119821) not detected in X-rays. The southern com- 
ponent is_a_binary itself with a projected separation of ~ 
0".13 dKohlerl l2008h . We have detected the [Ne II] line to- 
wards the T Tau system, co nsistent with previo us studies (e.g., 
Ivan den Ancker et all 119991: iGudel et al.l (20101) . Based on the 
analysis of the continuum-subtracted line, there is a contribu- 
tion from extended emission, in agreement with previous results. 
The first detections of [ Ne II], [Fe II], and [O I], am ong other 
gas tracers, come from Ivan den Ancker et alJ (|1999), based on 
Infrared Space Observatory (ISO) observations. However, we 
do not detect [Fe II] emission, which would be expected to be 
present in an outflow environment. Van Boekel et al. (2009) stud- 
ied this system (R~30 000) and have determined that the emis- 
sion is extended and associated with an outflow from the North 
and South components. Only a small fraction of the emission is 
related to the X-ray bright T Tau N. 

XZ Tau is a close binary system with separation o f 0."3. Both 
stars have collimated, bipolar jets (Krist et aD l2008l) . Due to the 
close separation, we cannot resolve the system with Spitzer. 
Observations have been made in 3 sets of 2 nod pointings; we 
detect [Ne II] only at the pointing with the highest flux in the 
continuum, centered at the position of the stars. This is the first 
detection of [Ne II] reported toward this source. 
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CoKu Tau 1 is a Class II obj ect. It is a binary with separation of 
0".24 dRobitaille et al.ll2007l) in which we have detected [Ne II] 
emiss ion. Previous observation s in Ha show evidence of an out- 
flow (Eisloffel & Mundtl[T998h . but we do not detect the [Fe II] 
emission usually observed in outflows. 

CoKu Tau 4 is a Class II object kno wn to have a large excess 
at ~ 30 pm and no excess at < 8 /mi (For rest et al.l l2004). This 
was interpreted as a clear signature of a tra nsitional disk, but re- 
cent observations dlreland & Krau s 2008) have shown that the 
clearing of the disk at shorter wavelengths is actually due to the 
presence of a stellar-mass companion. The binary projected sep- 
aration is calculated to be ~ 7.8 AU. We have detetected [Ne II] 
emission towards this binary. 

GM Aur is a Class II so urce with [Ne II] e mission; this object 
is also a transitional disk (ICalvet et al.1 2005). We have dedicated 
background observations, so we can discard the possibility of 
the emission coming from the background or a nearby source. 
The first detection of [Ne II] in the Spitzer spectrum was re- 
ported by ICarr & N aiita (2008). Foll ow-up observations of the 
[Ne II] line at R~80 000 dNaiita et al.ll2009b have confirmed that 
the em ission is from a disk rather than an outflow. Espaillat et al. 
d2007l) have derived a 5<x upper limit for the luminosity of the 
[Ne II] line of 7.4 x 10 28 erg s _1 . This result is based on ob- 
servations made with Spitzer IRS low-resolution modules (R~ 
60-120). Our calculations based on the high resolution module 

give a luminosity of (1.5 ± 0.3) x 10 28 erg s" 1 . 

RW Aur is an optically bright T Tauri star; [Petrov et"al] d2001l) 
suggested that this is a binary system. The secondary, RW Aur B, 
is located 1".2 away from the primary. We have detected both 
[Fe II] lines (17.93 and 25.99 /mi) in the spectrum of this sys- 
tem, which is not spatially resolved by Spitzer. The observations 
were made in 3 sets of 2 nod pointings. The line at 17.93 pm is 
detected in one of the offset positions, indicating that the emis- 
sion is likely to be extended. The [Fe II] line at 25.99 pm is 
detected in all the three positions: centered on the position of 
the star, and two in positions displaced from the star position. 
The contribution from extended emission is important, but the 
line flux is higher at the position of the source. Previous obser- 
vations have shown emission of [Fe II] in the near-IR from t he 
jet associated with this system dHartigan & Hillenbra nd 2009). 
UY Aur (A, B) is a binary system with separation 0".88 
dHirthet al.lll997h . Observations were made in 3 sets of 3 nod 
pointings. We have detected [Ne II] emission at the 3cr limit. 
The line is detected only in the pointing centered in the position 
of the system. We note that in this case the line is broad, with a 
FWHM 1.6 times higher than the IRS instrumental FWHM. This 
system is know n to drive a high v elocity (~ 300 km s~') [O I] 
bipolar outflow dHirth et alj|T997), We do not observe any shift 
of the position of the line; the broadening could indicate either a 
strong wind or a jet with a blue and red-shifted co mponent. 
MHO 1/2 is a binary with a separation of 3" (B riceno et alJ 
119981) . Since the Spitzer observations do not resolve the system, 
we cannot know in which star the emission originates. The sys- 
tem has been observed in 3 sets of 2 nod pointings. The [Ne II] 
line is detected only at one of the three positions. 
V773 Tau is a multiple system. V773 Tau AB is a WTTS spec- 
troscopic binary with a separation of 0".003, i.e.,~ 0.37 AU at 
the distance of Taurus dWeltvi l995). V773 Tau C is a visual com- 
panion at 0".17 from V773 Tau AB, classified as a CTTS. A 
fourth member of this system, V773 Tau D, has a n angular sep- 
aration of ~ 0".22 with respect to V773 Tau A dWoitasll2003l) 
and it has been cla ssified as an infrare d companion due to the 
shape of its SED dDuchene et alj|2003h . These objects are be- 
lieved to be T Tauri stars embedded in dusty envelopes and un- 



dergoing episodes of intense accretion. But their true nature is 
not yet understood. The Spitzer observations do not resolve the 
c omponents of the syste m, but according to the SED presented 
in iDuchene et all d2003h we assume that V773 Tau D is respon- 
sible for the bulk of the mid-IR emission and therefore for the 
level of continuum observed. However, we cannot be sure that 
the [Ne II] emission comes from it. In fact, it could also be pro- 
duced in shocks from one of the components or in the interstellar 
environment in a similar way as the case of T Tau. The stellar 
properties for the system are compiled in TableQ] The Ha emis- 
sion and X-ray luminosity probably come from V773 Tau AB, 
the WTTS binary, while the mass accretion rate value adopted 
is the one derived for V773 Tau C, which is likely similar to 
V773 Tau D. We assume that in the infrared regime, the emis- 
sion is dominated by the CTTS, so we have included the system 
in the group of Class II sources. In the correlation tests presented 
in Section [5] we have not used the Ha equivalent width, which 
is derived for the WTTS binary. 

HBC 366 is a WTTS. We have detected one of the H 2 lines 
(17.03 pm) in this spectrum. No other gas lines have been de- 
tected. 

IRAS 04239+2436 is a binary of 0".3 separation, driv ing a 
highly collimated near-IR [Fe II] jet dReipurth et alJ l2000). The 
observations were done in three sets of two nod pointings. We 
detect both iron lines (17.93 and 25.99 pm) in the IRS spec- 
trum. The two lines are detected in all three sets of pointings, 
but when analyzing the continuum-subtracted line, the highest 
contribution comes from the position centered on the star. The 
luminosity of the line at 25.99 pm is more than a factor of 3 
high er than the luminosity at 17.93 /mi. Previous observations 
dDavis et al.ll2001blDavis et al.ll2003l) have detected emission of 
H2 (2.122 pm) and [Fe II] (1.644 pm). The Spitzer spectrum of 
this source does not show any evidence of H2 emission. 
IRAS 04303+2240 is a CTTS. We have detected both [Fe II] 
lines (at 17.93 and 25.99 pm), with the line at 17.93 pm detected 
at the 3<x level. Although there is an excess at the wavelength of 
the [Ne II] line, it is formally not detected (at the 3<x detection 
limit). The ratio between the luminosity of both [Fe II] lines is 
1 :2. Observations were done in 3 sets of 2 nod pointings; the line 
at 17.93 pm is detected only in the position centered at the star. 
The line at 25.99 pm is detected in two out of three positions; the 
contributions of the centered position and the offset position are 
similar. Both lines are broadened by a factor of two compared to 
the intrinsic line width of the instrument. There is no information 
in the literature about previous detections of [Fe II] towards this 
source or about jets or outflows that could explain the [Fe II] 
emission, or winds explaining the broadening of the lines. 
V710 Tau is a binary system formed by a CTTS and a WTTS. 
Obser vations in the X-ra ys made with Chand ra ([Shukla et al. 
2008) and in the near-IR ( White & Ghezi l2001) have shown that 
the separation of the binary is 3". 2. We detect H2 at 17.03 and 
28.22 pm. Considering that the star responsible for the contin- 
uum at mid-IR wavelengths is the CTTS, V710 Tau A, we have 
used the Ha equivalent width and X-ray luminosity of that star in 
the present study, and included the system in the Class II group, 
but the lines could also be produced by the WTTS. 
SST 041412.2+280837 was included in a sample of the low- 
mass star population of Taurus for which we have obtained 
Spitzer-IRS spectroscopy. It was previously kno wn as IRAS 
0411 1+2800G and is classified as a Class I object dRebull et al.l 
120101) . This is the only object in our sample for which we have 
used low-resolution spectra (SL and LL) since we do not have 
the SH and LH spectra. We have detected [Ne II] and [Fe II] at 
17.9 and 25.9 pm towards this object. 
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SST 042936.0+243555 is a newly discovered member of Taurus 



Class II. The spectrum shows [Ne II] emission. 
SST 043905.2+233745 al so belongs to the sa mple of newly dis- 
covered Taurus me mbers dRebull et alJl2010l) . It is a Class I ob- 
ject as classifi ed by iLuhman et al.l (1201 Oh . but a Class II object 
according to Rebu ll et al.l d2010 !). We have detected emission 
lines from [Ne II] and H2 at 17.0 fim. 

Appendix B: [Ne II] non-detections. 
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Fig.A.3. Archival SH spectra of DK Tau and BP Tau showing a variety of emission lines fro m water and organic molecules. 
We present for comparison the spectrum of AA Tau prev iously presented in lCarr & N aiita (2008). The empty diamonds show the 
position of some water lines (see Pontoppid an et al.l2010h . The spectra were background subtracted and shifted to allow comparison. 



Baldovin-Saavedra et al.: Searching for gas lines in the young population of the TMC, Online Material p 9 



BP Tau 



CW Tau 



„ 0.47 

x 0.46 
0.45 
0.44 



0.42 
0.41 



12.6 12.8 13.0 

DN Tau 

0.36 [ } 
0.35 i 

0.34 ' ^ 

0.33 "VvjV^ 
0.32 i 

12.6 12.8 13.0 



GG Tau 



„ 0.66 
j>> 0.64 
" 0.62 
x 0.60 
2 0.58 
fc 0.56 



12.6 12.8 13.0 



HP Tau 




12.8 13.0 



UZ Tau 



1.10 
1.05 
1.00 
0.95 



12.6 12.8 13.0 

HD 35187 

4.5 I i 1 

4.4 
4.3 

4.2 

4.1 Vij^ 



1.45 
1.40 
1.35 
1.30 
1.25 
1.20 



12.6 12.8 13.0 
DO Tau 




0.100 
0.095 



12.6 12.8 13.0 
HD 36112 



3.2 
3.1 
3.0 
2.9 
2.8 



CZ Tau (AB) 
0.85 
0.80 
0.75 



12.6 12.8 13.0 

DR Tau 



2.2 
2.1 
2.0 



12.6 12.8 13.0 



GK Tau 



0.86 [ 

0.84 

0.82 



0.74 1 

12.6 12.8 13.0 



RW Aur A 



1.35 
1.30 
1.25 
1.20 
1.15 



12.6 12.8 13.0 



Coku Tau 3 



0.40 
0.39 
0.38 
0.37 
0.36 
0.35 



12.6 12.8 13.0 



Haro 6-28 



0.18 
0.16 
0.14 
0.12 
0.10 




DF Tau (AB) 

o.9o ; 

0.85 
0.80 
0.75 
0.70 



DK Tau 



DL Tau 




12.6 12.8 13.0 

FV Tau (AB) 
1.80 r 

1.75 
1.70 
1.65 
1.60 
1.55 



12.6 12.8 13.0 



0.31 
0.30 
0.29 
0.28 



FX Tau 



12.6 12.8 13.0 
GV Tau (AB) 

25 



12.6 12.8 13.0 



HL Tau 




RY Tau 



">.% 

7.4 
7.2 
7.0 
6.8 
6.6 

12.6 12.8 13.0 




12.6 12.8 13.0 



Haro 6-37 



13.0 



SU Aur 

2.8 
2.7 
2.6 • 
2.5 

12.6 12.8 13.0 



HBC 366 



078 
076 
074 
072 
070 
068 
066 



12.6 12.8 13.0 



V710 Tau 




0.90 
0.85 
0.80 



12.6 12.8 13.0 
FZ Tau 

0.95 . 
0.90 ] 
0.85 
0.80 
0.75 . 

12.6 12.8 13.0 



HN Tau 



0.74 
0.72 
0.70 
0.68 
0.66 
0.64 



12.6 12.8 13.0 

UX Tau A 




HD 31648 




12.6 12.8 13.0 



12.6 12.8 13.0 



12.6 12.8 13.0 



12.6 12.8 13.0 



12.6 12.8 13.0 



12.6 12.8 13.0 
V892 Tau 

30 



28 
27 
26 

12.6 12.8 13.0 



IC2087IR 



5.4 
5.2 
5.0 
4.8 



12.6 12.8 13.0 



04016+2610 



4.2 
4.0 
3.8 
3.6 



12.6 12.8 13.0 



F04101+3103 

2.3 
2.2 
2.1 



12.6 12.8 13.0 



04108+2803B 




12.6 12.8 13.0 



04239+2436 

2.8 I ] I 
2.7 i 
2.6 
2.5 
2.4 i 

2.3 1 j : 

12.6 12.8 13.0 



04278+2253 

9.0 



12.6 12.8 13.0 



F04299+2915 

0.18 




12.8 13.0 
A [/im] 



04303+2240 

2.0 r 

1.9 
1.8 
1.7 

12.6 12.8 13.0 



04361+2547 




04365+2535 



12.6 12.8 13.0 
A [/im] 




12.6 12.8 13.0 
A [^m] 



Fig. B.l. Spectra around the [Ne II] line of the targets in which we did not detect the line. The position of the line is plotted 
dashed-dotted line and the 3<x noise level is plotted in dashed line. The name of each star is labeled in the top of each plot. 



